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The Drell-Yan mechanism for the production of lepton pairs is one of the most basic 
processes for physics studies at hadron colliders. It is therefore important to have accurate 
theoretical predictions. In this work we compute the two-loop virtual mixed QCDxQED 
corrections to Drell-Yan production. We evaluate the Feynman diagrams by decomposing 
the amplitudes into a set of known master integrals and their coefficients, which allows us to 
derive an analytical result. We also perform a detailed study of the ultraviolet and infrared 
structure of the two-loop amplitude and the corresponding poles in e. 



I. INTRODUCTION 

The Drell-Yan process is one of the most precise probes available at hadron colliders. It allows 
for precise measurements of the gauge boson masses IMM, widths [|4J, ISfl and asymmetries [|6|] 
and is very sensitive to physics beyond the Standard Model like new gauge bosons jvl- lioll . One 
reason it is such a powerful probe is its very simplicity. Its experimental signature, two leptons 
plus anything, is quite robust against radiative emission. Theoretically, it is perhaps the simplest 
process to compute at hadron colliders. It was the first hadronic scattering process to be computed 



at next-to-next-to-leading (NNLO) in QCD |llll,ll2|], almost twenty years ago. 



In recent years, it has become clear that electroweak corrections lll3Ml6ll to Drell-Yan produc- 
tion are also very important. Electroweak corrections can distort the line-shape and thereby affect 
the measurement of the gauge boson masses. Radiative corrections can also become very large at 
the high energies (several hundred GeV) which will be probed at the LHC. 

An important next step to refining the prediction for Drell-Yan production is the calculation of 
the complete mixed QCD and electroweak corrections. Currently, only the virtual corrections to 
the quark - gauge boson vertex are known in the literature [17]. We embark on this project by 
computing the simplest gauge-invariant part, the mixed QCDxQED virtual corrections. That is, 
we ignore all W and Z boson interactions, and consider only virtual photon and gluon exchanges. 
In addition, we take all fermions to be massless (except the top quark, which does not enter into 
this part of the calculation). For most of the calculation, there is no barrier to including a non- 
vanishing lepton mass. For the box contributions, however, one would need new master integrals 
with two massive external legs and massive internal propagators. The effects of non-vanishing 
masses (at least for components that do not involve box contributions) will be added at a later 
stage of the project. 

The outline of this paper is as follows. In Section |ll] we define some generalities and our 
notation. In Section |lll]we give an outline of the calculation and in Section|IV]discuss the structure 
of the ultraviolet and infrared poles. The results are presented in Section |V] and we present our 
conclusions in Section |Vll The Appendix contains supplementary information about the next-to- 
leading order process as well as about the master integrals arising in the calculation. 



II. GENERALITIES AND NOTATION 

We study the Drell-Yan process of quark (q) anti-quark (q) annihilation into a charged lepton (£) 
pair 

qipi)+q{p2)^riP3)+i+{p4) (1) 

where Pi^PiiPSiPA denote the momenta of the particles, which are all considered as incoming 
with Pi+ P2+ P3+ P4 = 0. In the following we will use the Mandelstam variables 

s = {Pl+P2)^ = {P3+P4f, 
t = {pi+P3)^={p2+P4f: 
U = (pl+p4)2=(p2+P3)^ (2) 

or Sij — {pi -\-pjY io express scalar products of the external momenta. 
The differential cross section is given by 
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color 



where the symbol A^^ denotes the number of colors of SU(A^c) and ^ is the matrix element. 
Within this work we consider only virtual corrections to the Drell-Yan process. The perturbative 
expansion of the corresponding squared matrix element is given by 
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color 



where Qq and Qi are the electric charges of the initial state quarks and the final state leptons 
in units of the elementary charge e. The symbols a and a, are the fine structure constant and 
the strong coupling constant, respectively; Cp = {N^ — \)/{2Nc) denotes the Casimir operator 
of the fundamental representation of SU(A^c)- The dots stand for higher order corrections. Here 
and in the following we will write the expansion of any function of a and oty as /(a, a,) = 
Lm,n{0i/n)"' ( a,/ :;r)" /("'"). The well-known leading order resuk A^^'O) of Eq. @ in J = 4 - 2e 
space-time dimensions reads 
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The one-loop QCD [| 181- 12 ID and one-loop QED corrections [|13|1 are known. For completeness we 
will give the bare results for A^^'*^) as well as for A*-*^'^) in Appendix lAl since they are needed for 
the subtraction of the ultraviolet poles through the renormalization procedure as well as for the 
identification of the infrared poles. 



m. CALCULATION OF THE BARE PROCESS 

Depending on the nature of the electromagnetic corrections, the QCD x QED corrections to 
Drell-Yan production can be broken up into four classes, which are: initial state corrections, final 
state corrections, mixed initial and final state corrections, and vacuum polarization corrections. 
Sample diagrams for each of these four classes are shown in FIG[il 
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FIG. 1 : Example diagrams which contribute to the two-loop mixed QCD x QED corrections. The spiral lines 
denote gluons, the wavy lines are photons and the straight lines are fermions which can either be quarks in 
the initial state or leptons in the final state. 



The initial state electromagnetic corrections consist of two-loop corrections to the quark - 
photon vertex. A sample diagram of this class of correction is shown in FIG. [IJa). All of the 



diagrams that appear in this portion of the calculation are topologically identical to diagrams that 
appear in the two-loop QCD corrections to Drell-Yan production. By simultaneously computing 
the two-loop QCD corrections and verifying the known result f22\, we obtain a strong check on 
this part of the calculation. We also include the interference of the one-loop QCD correction and 
the one-loop QED initial state correction. There is but one diagram of each sort. 

The final state virtual corrections are quite trivial, since the only contributing two-loop diagram 
is the one shown in FIG. \V[h), which is just the product of two one-loop triangle diagrams. The 
virtual corrections to this channel also get a contribution from the interference of the one-loop 
QCD corrections and the one-loop final state QED corrections. Again, there is but one diagram of 
each sort. 

The mixed initial and final state electromagnetic corrections are the most complicated terms in 
this calculation and the only ones which involve the kinematic variables t and u in the loop inte- 
grals. A sample diagram is shown in FIG. [He). Even the interference of the one-loop amplitudes is 
relatively complicated as the electromagnetic part involves the sum of two one-loop box integrals. 

The vacuum polarization correction terms are also easy to compute as the loops are simple two- 
loop propagator integrals, like that shown in FIG. [Hdi), or the product of a one-loop propagator 
integral and a one-loop triangle as in FIG.[Hd2). The interference of one-loop amplitudes is again 
very simple as it only involves a single vacuum polarization diagram interfered with the single 
one-loop QCD diagram. 

We have performed two independent calculations of the virtual corrections and find complete 
agreement. The Feynman diagrams are generated with QGRAF [23]. The symbolic algebra pro- 



gram FORM [|24l1 is used to implement the Feynman rules, to interfere the two-loop diagrams 
with the tree-level contribution and to reduce the result to a set of Feynman integrals to be deter- 
mined. The calculation proceeds in two steps. In the first step all the loop integrals are mapped 
onto a small set of master integrals with the traditional integration-by-parts (IBP) method [|25I1 



in combination with Laporta's algorithm |[26L[27[]. In the second step these master integrals are 
evaluated. 

In one calculation, the integrals are reduced to master integrals using the program RE- 
DUZE [|28n . In the second calculation, the reduction has been performed with a FORM [|24ll29l . l30l] 



based implementation which uses the packages Q2E and EXP [3 ll 13211 to identify the different 
topologies and to adopt the proper notation. The program FERMAT [|33ll is used to simplify the 
rational functions in the space time dimension d. 

We find that at the two-loop level, all integrals can be expressed in terms of eight master topolo- 
gies which are shown in FIG. [2l All of the needed master integrals are known analytically in the 
literature to sufficiently deep order in the £ expansion and will be discussed in more detail in Ap- 
pendix [Bl At the end, the reductions and the evaluations of the master integrals are substituted 
back into the FORM program to produce the final result. 

IV. THE POLE STRUCTURE OF THE PROCESS 

A. Ultraviolet Structure and Renormalization 

We have performed renormalization in the MS scheme. Since we treat all particles as being 
massless, we only need to renormalize the couplings. Coupling constant renormalization is gov- 
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FIG. 2: The master topologies of the two-loop calculation ordered according to the number of internal lines. 
The arrow on the lines denotes the momentum flow. The symbols ki and ^2 are loop momenta. The integral 
representation of diagram (^2) has an additional irreducible scalar product (^1 + ps)^ in the numerator 
indicated by the (g) symbol. 
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where /i is the renormalization scale, Nu is the number of up-type quarks, A^^ is the number of 
down-type quarks and Ni is the number of charged leptons, while Qu, Qd and Qi are their electric 

9 1 

charges, +f , — ^ and —1, respectively. The symbol Ca = Nc denotes the Casimir operator of the 
adjoint representation of SU(A^c) and 7/ = 1/2 is the normalization of the QCD charge of the 
fundamental representation. The bare and renormalized couplings are related by 
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where e ~ 2.71828 is Euler's number and Ye — 0.577216 is the Euler-Mascheroni constant. Since 
the leading-order contribution to the squared matrix element is of order a^ and we compute 
terms through order a^ a,, we need to keep QED renormalization terms proportional to j8 ^^'^^ 
and /3^^'^^ [l34ll . while QCD renormalization does not contribute to our result. The fine- structure 
constant can be converted, if needed, from the above MS scheme to the on-shell definition with a 



conversion factor. This conversion is known to four-loop order in QED [|35L 13611 



B. Infrared Structure 

An important check on our calculation is to verify that we have obtained the correct infrared 
structure. Some years ago, Catani [37] proposed a formula predicting the leading poles (£~^ 
through e^) of two-loop QCD amplitudes. At that time, the e^ poles were presumed to be 
process dependent and therefore unpredictable. Nonetheless, direct calculations [38-42] showed 
that the e^ terms seemed to follow a simple pattern based upon the numbers of quarks and gluons 
that made up the external legs of the amplitude. 

Subsequently, Sterman and Tejeda-Yeomans [43] reformulated Catani's observation and identi- 
fied the origins of the various terms. They also identified the then-unknown term, the second-order 
correction to the so-called "soft anomalous dimension" which prevented the prediction of the e~^ 
terms. Aybat, Dixon and Sterman ll44l 14511 have since computed the two-loop corrections to the 
soft anomalous dimension, permitting the prediction of the full infrared structure of two-loop QCD 
amplitudes. 



C. The Infrared Structure of QCD Amplitudes 



For a general 2 — )• n scattering process, 

/l (Pl, Cl) 4-/2(^2, C2) -> /3(;?3, C3) H \-fn+2iPn+2: Cn+l) 



(10) 



where /, represent the flavors of the partons, pi their momenta and c,- their colors, we can write the 
amplitude as a vector in the space of color tensors {(Q) r^..}} as BVL kg. 14711 
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- e) ) = £^f,L (^pu ^, a,(ju2), e j x (Cl){,.} . 



(11) 



where 2 is an (arbitrary) overall scale and /i is the renormalization scale. 

In the formulation of Refs. 1.43- 45.1 . a renormalized amplitude may be factorized into three 
functions: the jet function ^f, which describes the coUinear dynamics of the external partons 
that participate in the collision; the soft function Sf, which describes soft exchanges between 
the external partons; and the hard-scattering function \Hf), which describes the short-distance 
scattering process 
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The notation indicates that |//f) is a vector and Sf is a matrix in color space. As with any factor- 
ization, there is considerable freedom to move terms about from one function to the others. It is 
convenient [441 1450 to define the jet and soft functions, ^f and Sf, so that they contain all of the 
infrared poles but only contain infrared poles, while all infrared finite terms are absorbed into |//f). 



1. The Jet Function 



The jet function ^f is found to be the product of individual jet functions ^f. for each of the 
external partons, 

^f(a,(Ai2),e) =n ^^■(a.(M'),e) • (13) 



Each individual jet function is naturally defined in terms of the Sudakov form factor [431], 
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The all-orders expression for the square root of the Sudakov form factor is 1148145 ill 
7/(a,(/i2),£)=exp<^-y -^ Xi{as{ii^),£)+'^i(-l,ds(^^,as{^'^),e^,e 
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The functions ^, ^ and yxi are anomalous dimensions that can be determined from fixed-order 
calculations of the Sudakov form factors for quarks and gluons [|22l l52l457n. Note that yia is the 
cusp anomalous dimension and J(fj is determined, order by order, from jKi- While the ^ are pure 
pole terms, the ^/ contain terms at higher order in e. 

The jet functions ^f. keep only the infrared poles from the logarithm of the form factor. The 
expansion of the jet function to second order in oc^ is 
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Nq is the number of quark flavors and ^„ = E^^j 1/fc" represents the Riemann zeta-f unction of 
integer argument n. The coefficients of the /3 -functions are given in Eqs. ([6117]). Even though the ^■ 
have terms at higher order in e, we only keep terms in the expansion that contribute poles to In ^,. 

2. The Soft Function 

Like the jet function, the soft function can be defined in terms of eikonal amplitudes and is 
determined entirely by the soft anomalous dimension matrix r5,, 
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In the color-space notation of Refs. [I37l. l46l l47h. the soft anomalous dimension is given by [|44l . l45h 
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where K = Ca (67/18 — ^2) — ^OTfNq/9 is the same constant that relates the one- and two-loop 
cusp anomalous dimensions. The T, are the color generators in the representation of parton i, 
multiplied by ±1, depending on the whether the parton is a particle or antiparticle and whether 
it is incoming or outgoing. In particular, outgoing quarks and gluons and incoming anti-quarks 
are multiplied by +1, while incoming quarks and gluons and outgoing anti-quarks are multiplied 
by —1. The conservation of color-charge is enforced by the identity Y,i T, = 0. Another useful 
identity is that T, • T,- = Q. 



D. The Infrared Structure of QED Amplitudes 

It has been found that the same factorization described in Eq. (fT2l) can be applied to pure QED 
amplitudes [421 1580. The two-loop amplitudes for Bhabha scattering and for e+e^ — )■ 77 in mass- 
less QED were found to obey the factorization formula of Catani [37] once the proper adjustments 
are made to transform the QCD anomalous dimensions into QED anomalous dimensions. 



The changes are as follows. The factors of the adjoint representation Casimir, Ca, originate 
from the gluon self interactions. As photons have no self interactions, Ca is set to zero. The 
fundamental representation Casimir, Cf is replaced by the squared electric charge of the fermion, 
Cf ^ QJ. The factors of TfNq originate from inserting fermion bubbles into the gluon propagators. 
In QED, the different types of fermions would be weighted by the squares of their electric charges, 
TfNq — 7- Nc-Nu Qu +NcNd Q^ +Ne Qj. In the soft anomalous dimension, the color charge matrices 
T, are replaced by the (scalar) electric charges Qj. With these changes, the anomalous dimensions 
for the QED jet function are 
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while the QED contribution to the soft anomalous dimension is 



(20) 



Using these parameters, one can predict the infrared structure of two-loop QED amplitudes, 
where the analog rules as in Eq. (fT9l) apply for the signs. When comparing to the results of 



Refs. 1421 15811 . one must account for the fact that those calculations are in the context of pure 
QED, involving only leptons and photons. As the universality of the £^^ terms had not yet been 
established, the (color diagonal) H^-^> factors for electrons and photons was quoted as 
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where A'^^ =NiQj. Transforming the results above into the notation of Ref. BTN . we find that the 
//(^) terms may be more generally written as 

,(2)_ ^4/^3 .. ^.. V ^2ft(2,0)/25 3 



w(2)_^/o(2,0)\2_ (3,0) 

where the subscript / indicates any charged fermion - lepton or quark. With these modifications, 
we find complete agreement with the results of Refs. [|42l. l58|l . 



E. The Infrared Structure of QCD x QED Amplitudes 

The leading terms in the infrared structure of QCD x QED corrections will come from the 
overlap of the one-loop terms for pure QCD and pure QED. The intrinsically QCDx QED terms 
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will be second-order contributions to the jet and soft functions. Based upon the way the parameters 
were determined for QED, we can make conjectures about the parameters for QCDx QED. Since 
the generators for photons and gluons commute, we should again set the Ca terms to zero. We 
need to be a little more careful about the Nf terms, however. Our approach is to tie the A^^ terms 
to the coefficients of the /3 -functions. The reason for this is that when the Nf term is part of the 
leading term in a /3 -function, it represents the insertion of a fermion bubble into a gauge boson 
propagator. Because the charge matrix of QCD is traceless, the bubble cannot connect a photon 
to a gluon and therefore these terms cannot contribute to a second-order mixed correction. When 
the Nf term is part of a second-order term in a j8 -function, however, it represents a term like 
those shown in FIG. [H which can represent a second-order mixed correction. Examining the two- 
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FIG. 3: Mixed second order contributions to the QED and QCD j3 -functions. 



loop anomalous dimensions in Eqs. (flTI) and (fT9l) , we see that the second order corrections to 
the cusp and soft anomalous dimensions are proportional to ^ = C^ (67/18 — ^2) — ^OTfNq/9 ~ 

(2/3 — ^2) Ca + 10/3 Pq^j^. Since we have argued that neither non-Abelian nor first-order /3- 
function corrections can contribute to second-order mixed corrections, we conclude that there are 
no mixed corrections to the cusp and soft anomalous dimensions at this order. That leaves only 
the ^- terms. By the same reasoning as for the cusp and soft anomalous dimensions, we set the Ca 

and A'^ terms to zero in forming ^^ ' , but we predict that the Cp term should be transformed into 

Cf QI- For ^^\y , we again drop the non-Abelian and first-order /3 -functions, but we predict that 

2fi(^^2)and^(^'^)-2fi(2'^) 



we should keep the second-order j8 -function terms to obtain ^i ' 



We can thus write combined expressions for the jet and soft functions which we claim are valid 
through second order in both QCD and QED, 
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F. The Infrared Structure of the Drell-Yan Amplitude 



We can now examine our result for the Drell-Yan amplitude to see if we match the expected 
infrared structure. We start from the factorization formula, Eq. (fT2l) . and expand both sides in 
powers of a and a,, 
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Because of the trivial color structure of the Drell-Yan amplitude, the soft anomalous dimen- 
sion matrix is proportional to the unit matrix and may be treated as a scalar function. The 
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squared matrix element of Eq. dH) is related to the decomposition of the amplitude in Eq. (|26]) 
by Isp.n I^P = {e^E/{An)f^ (f )^ {^dy \J^dy)- The values of the jet and soft functions for the 

color 

Drell- Yan process are given by 



^i"'-(2^-^)te^e?)' ^^"-(2^-^) 



DY 



1 



+ ^^ + 



4e4 4e3 16£2 



c4e5^e,2)_^P _3^^^3^^|c^g2 



c(l,0) _ 1 

^^^ ~ Ye 



(ej+eain(^)+2e,e,(in(^)-in(^ 



;(0.1) 

'or 



^DY 



2e V-5 

4£2 V -5 



(ej + e?) In ( M +22, e^ (in ("^ 1 -in 



M^ 



(28) 



We find complete agreement between our result and the expected infrared structure presented in 
Eq. dlH]), including the intrinsically QCDxQED term in ^jyy^^ ■ 

V. RESULTS 



As our final result we present the interference of the finite hard-scattering terms that appear in 
Eq. (127]), defined by 



2('^)'Re 

.nJ 



H 



(1,0) 
DY 



//XVn + (//^^^ 



^DY 



^DY 



H 



(2,0) 
DY 



N,QlQJe''CFB^'''^ 



(29) 



where we performed the renormalization in the MS scheme as described in Section IIVAI The 
infrared poles are subtracted in d dimensions with the help of Eqs. (I27H28I) . We decompose this 
mixed QCD x QED two-loop contribution with respect to the charge factors 



B 



"•"^e,&s<i'"+^ 






(30) 



where the sum over i' and q' runs over all leptons and quark flavors which are active in the closed 
fermion loop. Each of the five terms corresponds to one of the classes of diagrams shown in 
FIG.[I](a)-(d) and corresponds to a gauge invariant subset of diagrams in this decomposition. The 
individual terms of Eq. (|30l) are 






(31) 



B 



(1,1) 



128-H^;r2 + ^;rV(l4;r2-96)log(^) 
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. |50-^.^)log^(i,)-.21og3(-l).2V(-L 



^(11) 155 140 2 ... /^28 2 92 , , 



's?' 



(^) 



112, 2/^* ^ 8, 3 



M' 



(1,1) _ 320 140 2, /28 2 104 N/ . 



112 



5 



(1,1) 



4i:Ij1lu{ — \-4 



M 



- + 4-log — 

5 5^ V 5' 



2,+ 9 



log' 



f«' 



s 



3a2 + ?2 



log 



s 
3 ^^^ -2- log ( — ) -2- log ( - 



s^ \ s 

f2 I 'if,, I i/;,,2 



19r + 3?M+16M^ \ , ( -u 
log — 



+ 21ogM:^ 



log^ 

-3 



+ 2 log, 2 



s 

t-u ^ .. 

log 

s \ s 



M^ 



log^ — +2 . log — 






5 
— M 



(32) 



(33) 



7 t^ + U^ 2 
^i — TT 



2 log 



—u 



+ %' 



4 ?-M 8m+15? ( -u 

35 35 \5/6 5^ 



5 m2 + 3^2 



log^ — 

5 



A 1 t — U 

TV --■ 



15 5 



3m^ + ?^, a f —u\ 2t + 2u, -x I —u 

— -r- i°g I — ) - o — log I — 



652 



5 y 3 5 



4 log2 — 

5 V 5 



3<-M,„,Ml|_„^„) 



(35) 



where Li„(z) = Y^=\ fn is the poly logarithm function and the symbol {t ^ u) stands for the same 
terms a given before only with the Mandelstam variables u and t interchanged. As an additional 
check of our calculation we have kept the complete dependence of the gauge parameter in the 
gauge boson propagators and have verified their cancellation. 

VI. SUMMARY & CONCLUSIONS 



We have computed the two-loop virtual corrections to Drell-Yan production at order asO? . 
The calculation of these mixed QCD x QED corrections includes two-loop corrections to the quark 
vertex, one-loop corrections to the quark and lepton vertices, vacuum polarization corrections 
to the photon propagator as well as two-loop box diagrams connecting the hadronic and leptonic 
states. The computation is accomplished by reducing all Feynman integrals to a small set of master 
integrals. The latter ones are known analytically to sufficiently high order in the £ expansion to 
allow us to derive an analytical result for the finite amplitude. 

We have also shown that the infrared structure of the mixed amplitudes follows from the same 
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universal factorization structure that governs the pure QCD and QED amplitudes and we have 
determined the value of the two-loop mixed anomalous dimension. 

Acknowledgments: We would like to thank Andreas Scharf and Doreen Wackeroth for useful 
discussions. This research was partially supported by the U. S. Department of Energy under 
Contract No. DE-AC02-98CH10886. 



Appendix A: Bare next-to-leading order results in terms of master integrals 

We present our bare results for the next-to-leading order processes in terms of the master inte- 
grals and coefficients to all orders in e. For the one-loop QCD and QED corrections we adopt the 
decomposition of the squared matrix element as given in Eq. dU; all quantities are considered as 
bare. We find 



,(0.1) 



a(O'O) / 2 
^fl---2e)55W: 



(Al) 



I (1,0) 



A),^- = (eJ+el)Aif'V £e^,^^..£e2 \a(o,o) 



(' 



l-£ 

2e-3 



B'oi^) 



+ QqQ( 



10---2e') ^-J^B'o{s) + (6-- + 2£) ^-^B\,{u) 



, u(t^ + 3u^) 1 ,■/ X / 

+ ( 2^ --3esu D'Q{s,u)-{t^u) 



(A2) 



with the integrals 

B'q{s) = {47t^^Y e-^^^2Re [^^^WJ and D'q{s,u) = {47t^^Y e-'=^^2Re \li^\s,u) 



(A3) 



where e ~ 2.71828 is again Euler's number and Ye ~ 0.577216 is the Euler-Mascheroni constant. 



.(1) 



The values of the master integrals I2 and I^ are given in Appendix IB In the coefficients of the 
master integrals of Eqs. (IA1HA2I) . spurious poles in e appear, which arise while solving the linear 
system of IBP equations. As a result, one must know the master integrals which are multiplied 
by such spurious poles at higher order in the e expansion. The same situation also occurs in the 
two-loop arnplitude. In principle, those spurious poles could be avoided by choosing an epsilon 
finite basis [|59n . However, since all necessary master integrals are known either in closed form or 
to sufficiently high order in e, we retain the standard basis of master integrals. 



Appendix B: Master Integrals 

The reduction process relates complicated integrals with many terms in the numerators and 
denominators to "simpler" integrals with fewer terms in both numerators and denominators. In 
general, it is preferred that the master integrals have numerators equal to unity, and denominators 
which only contain propagators of unit strength, but this preference cannot always be satisfied. 

In this calculation, we encounter eighteen two-loop master integrals. Of these, eight represent 
distinct topologies which are shown in FIG. |2l the others are related to these eight by relabeling the 
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external legs. Only one of the distinct topologies has an irreducible numerator (or, equivalently, a 
doubled propagator in the denominator). 

All of the master integrals needed for this calculation are known in the literature. The double 
box integrals, FIG. [2lgi) and FIG. WLg,2) are known as Laurent expansions in the dimensional 
regularization parameter £. The others are all known in closed form and can be readily computed 
using standard Feynman parametrization techniques. 

(k) 

In the following we define the master integrals Ip-J with loop momenta k^ and ^2 in Minkowski 
space, where the superscript k indicates the number of loops, the subscript p denotes the number 
of propagators and s enumerates integrals with the same number of loops and propagators. For 
clarity, we also indicate the Mandelstam variables that appear as arguments. 



1. One-loop master integrals 

At one-loop order we have the five master integrals /j (s), /j (t), I2 (u), I^ {s,t) and 
/j (5, w), which are define by 

2^^-' Ji7td/^D,D,' U{^.u)-e J .^d/2D,D2D3D4' ^ ^ 

with 

Di=k\ + ie, D2 = {ki-pi)^ + i£, 

D3 = {ki-pi-p2)^ + i£, D4 = {ki-pi-p2-p3)^ + i£. 



Their values are 



,.,,-_, .,..rr_H^0^^ 



4V)^^Q^pt^%£l ,_.)-.,^,,,,_e;,_e;,^ 



St r(l-2e) 



s 



t 



+{-t)-^2Fi{l,-e;l-e;l + -^ 



(B3) 



where 2F\{a,b;c;z) = lLk=oi'^)k{b)k/{c)kZ^ /k\ are hypergeometric functions, (a)„ = T{a -\- 
n)/Y{a) is the Pochhammer symbol and T{x) is the gamma function. 

2. Two-loop master integrals 

a. Three-Line Topologies 
There is one distinct three-line topology, shown in Fig.[2ta), which we label l\ \s) and define 



lf\s)=e^^y^ 1^^^^^, (B4) 



by 

with 

D5 = {ki-p2)^ + ie , D(, = iki-k2)^ + i£ , Dt = {k2 + pif + ie. 
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Its value is 



2)/.^ „„^ ,(2), 



/n.) = .^"^(-.) — - Vg:3,v""' - <^^> 



In addition, we also need /^ (?) and /^ (m). 

fe Four-Line Topologies 

(2\ 12') (2) ^2) 

There are four four-line master integrals, I\.{{s), I\.2{s), ^4-2(0 ^^^ m-2(") ^^^^ ^'^^ distinct 
four-line topologies. One is a simple iterated bubble diagram, shown in Fig. [2b), which evaluates 
to the square of the expression in Eq. (|B2I) : the other is shown in Fig. Oc). They are defined by 

/(2)f,l =,2arE / ^'^l ^^"^2 1 (2). . ^ 2ay, / ^^^1 ^^"^2 1 

(B6) 
with 

D8=kl + i£, D9 = {k2-pi-p2)^ + i£, Dio = {ki+pi)^ + i£, 

and are given by /j.j (5) = (/2 (i'))^, 

/2).. _2ey,( 2a ^1 - 2£) FVe) r(l + g) r(l + 2£) 

^4;2^'^ - ^ ^ "^ 2(l-2e)r(2-3£) • ^ ^ 

c. Five-line Topologies 

(2) (2) (2) (2) 

There are six five-line master integrals I^.({s,t,u), I^.({u,s,t), I^.l{s,u,t), /^.{(f,^, m), 

(2\ (2) 

/j. 2 (*,?,«) and 1^.2{s,u,t) with two distinct five-line topologies, which are shown in Fig. [2d) 
and Fig. Wi^)- The first topology has a bubble connecting two adjacent comers of a box, the other 
five-line topology is a box diagram with a diagonal-line connecting opposite comers. They are 
defined by 

5^1^ ''^ " J iTf^l^ iTf^l^ D,D,DeDnDn ' ^ ^ 

d'^ki d''k2 1 

Dn = (h +p3)^ + i£, Dn = {h - P2f + i£ -. 



(2), , 9PT/ /■ d'^h d''k2 1 



with 

Their results read 



Iflis^uj) 



e 



2ey£r2(-e)r(-l+2e) 



ni-3e) 



(-M)-2£r(l-e)2Fi(l,-e;l-e; 



+ (-5)-2^r(l + e)r(l-2e)2Fi(l,£;l-£;-^) 



5/ 
(BIO) 
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,P)„,„,,) ,,2... r3(-s)r(2.) 



(-") 



-2£ 



2fr(l-3e) 
+ i-s)-^' (l -2F1 (1, -2e; 1 -2e; -^ 



I-2F1 l,-2e;l-2£; 



(Bll) 



d. Six-line Topologies 

There is only one six-line master integral, the non-planar triangle diagram, shown in Fig. [2f) 
and defined by 

lj^\s) =e^'r^ f^^ ^ , with Dn = {k2-h-p2f + ie. 

Its result can be expressed with the help of generalized hypergeometric functions 
pFq{au. ..,ap;bi,...,bq,z)= 1^=0(^1)^ • • • K)fc/((^i)fc- • • ipq)k)z''/k\ by 



/f(^) 



-e2^^£(-5)-2-2er(l + 2e) 
r(l-e)r4(l-2e)r(l+e)r2(l+2e) 4r2(l-£)r(l-2£) 



e4r2(l-4e)r(l+4e) 



e4r(l-4e) 



r2(i-£)r(i-2£)r(i + £) 

2e4r(l-3e) 

4r2(i-£)r(i-2£)r(i+£)r(i+2£) 
£4r(i-4£)r(i + 3£) 



5^2(1, -2£,-4£;l-2£,l-3£;i; 



3F2(£,£, l+2£; l+£, l+3£; 1) 



r3(i-£) 
■2£4r(i-3£; 



■4^^3(1, l-£, -2£,-4£;l-2£, l-2£, l-3£; 1; 



(B12) 



Note that the closed-form expression given above appears to differ slightly from that given by 
Ref . [60] . However, by rearranging the F-functions and applying various hypergeometric identi- 
ties, one finds that the two expressions are exactly equal. 



e. Seven-line Topologies 



ri^), 



t(^)i 



Ti^)l 



r(2)/ 



There are four seven-line master integrals Ij.({s,t), l).{{s,u), l).2{s,t) and I-j.^is^u) with two 
distinct topologies. One is the double-box topology where all propagators are of unit strength, 
shown in FlCOgj). It is defined by 



I).{{S,U) =6^'-^^ / T77 77= 

7,1V ' ^ J .a/2 ij^d/z 



1 



d/2 DiD3D5DeD8D9Di4 



with Di4 = {k2+P3f + ie, (B13) 



and known as a Laurent expansion in £ n6u\ 

(-5)-2-2e 



4?i(^,") 



^+44 [2^'-^^^^]-^ 



4Li3(-x)-4£Li2(-x) 
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+ 2Lii(-^)(£2 + 6C2) + ^^3 + 33C2^-yC3]+^^' + 36C2^'-yC3^ 

+ 87 ^4 -4 (S2,2 {-x)-iSi2 (-x)) +44Li4 (-;c) +4Li3 {-x) (Lii {-x) - 6 
+ 2U2{-x) {f-2iLk{~x) + 20i;2)+Ul{-x) (f + eCi) 

X 



(B14) 



with i = log(;c), X = u/s and the generalized polylogarithm function S„,p(z) = 

{-ir+P-^/(n-l)\/p\jidt'\og"-' {t')logPil -zt')/t'. 

There are two equivalent representations for the second seven-line topology. One is the double- 
box with a doubled propagator. The other representation is a double-box with an irreducible 
numerator, shown in FIG. |2tg2)- The latter is defined by 



I^lis^u) =e'''' I 



2eye fd%_d%_ 



ih+P3] 



in'^l^ in'^l^ DiD^D5D(,D^D>)Du ' 



(B15) 



When one uses the integral with the doubled propagator, the reduction procedure generates a spu- 
rious pole in £, meaning that one needs the double-box integrals expanded to order £^. When one 
instead uses the above double-box with an irreducible numerator, the reduction does not generate 
the extra pole, meaning that one only needs to expand the integrals to order e*^. 

The double-box with an irreducible numerator was first calculated in [|62D, with the result 



I^^l{s,u) = e^^^^ 



i-s) 



-2-2e T-2 



r^(i+e) 



4 £4 



14C2 , 1 



^£3 + 28^2^ 



+ 4{f + 6^2)Ui (-x) +8Li3 i-x) - 8£Li2 {-x) - 16 Cs] - -^"^ -26 Ci^^ 



■y^^ + 52C2 



Lii(-Jc)-5 [f + 6i;2] Lil i-x) + [6 f + 20 iUi{-x) 

- 8 C2] Li2 i-x) + [Si- 20Lii {-x)] Lig {-x) + 20 82,2 {-x) - 20i 81,2 {-x) 

- 28Li4(-x) + [28£-F20Lii(-x)] ;3-14C4 + ^(e) }■ (B16) 
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